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ABSTRACT 

Cool and dense prominences found in the solar atmosphere are known to be partially ionized because 
of their relative low temperature. In this Letter, we address the long-standing problem of how the 
neutral component of the plasma in prominences is supported against gravity. Using the multiple 
fluid approach we solve the time-dependent equations in two dimensions considering the frictional 
coupling between the neutral and ionized components of the magnetized plasma representative of 
a solar prominence embedded in a hot coronal environment. We demonstrate that given an initial 
density enhancement in the two fluids, representing the body of the prominence, the system is able to 
relax in the vicinity of magnetic dips to a stationary state in which both neutrals and ionized species 
are dynamically suspended above the photosphere. Two different coupling processes are considered 
in this study, collisions between ions and neutrals and charge exchange interactions. We find that for 
realistic conditions ions are essentially static while neutrals have a very small downflow velocity. The 
coupling between ions and neutrals is so strong at the prominence body that the behavior is similar 
to that of a single fluid with an effective density equal to the sum of the ion and neutral species. We 
also find that the charge exchange mechanism is about three times more efficient sustaining neutrals 
than elastic scattering of ions with neutrals. 

Subject headings: plasmas — magnetic fields — Sun: corona 


1. INTRODUCTION 

Traditionally the problem of solar prominence support 
is described in terms of the magnetic force that balances 
the solar gravitational force. In the existing models it 
is assumed that the prominence plasma is fully ionized. 
However, observations of the prominence body in Ha, 
which is an excitation line, suggest that the observed 
plasma can not be fully ionized. The same happens for 
prominence observations in the He I 10830 A line. At 
a temperature of 20, 000 K, typical of the prominence 
corona transition region (PCTR), we already have ion¬ 
ized H, ionized He and also ionized Ca. Thus, the sup¬ 
port of both the ionized and neutral components is re¬ 
quired. The main problem is that the neutral component 
is not directly affected by the restoring magnetic forces 
that can balance gravity. 

In the past it has been proposed that the frictional 
coupling between neutrals and ions through elastic col- 
lisions might play a role in the support of prominences. 
iMercier &: HevvaertsI (|1977ll examined the relative down¬ 
ward diffusion of neutral atoms due to gravity and tak¬ 
ing into account the magnetic field. These authors found 
that the resulting downward velocity is irrelevant to ex¬ 
plain the mass loss in prominences. iBakhareva et all 
(|1992f) did the first attempt to construct a model based 
on the one-dimensional Kippenhahn-Schliiter solutions 
including ion-neutral collisions and provided some de¬ 
tails ab out the dynamical regim es of prominence evo¬ 
lution. iPecseli fc Engvolil (I2n0(t l suggested that wave 
damping caused predominantly by ion-neutral collisions 
in the p rominence core may balance the acceleration of 
gravity. [Gilbert et al.l (l2002[ l have shown that the drain¬ 
ing effect for a hydrogen-helium plasma is rather small, 
specially for the hydrogen that moves down in the di- 

jaume.terradas@uib.es 


rection of the photosphere at a typical velocity of only 
3.7ms“^. The results of these works indicate that ion- 
neutral coupling may be quite relevant to sustain neu¬ 
trals, nevertheless, a consistent study taking into account 
frictional coupling, magnetic forces, and gravity has not 
been developed so far. 

Another mechanism that provides a frictional force is 
the resonant charge exchange process. Contrary to the 
elastic scattering between ions and neutrals, charge ex¬ 
change collisions are not identity-preserving. For hy¬ 
drogen an energetic proton captures the electron from 
a lower-energy neutral producing an energetic neutral 
which has essentially the same energy as the incident 
proton . The reader is referred to lGoldst on fc R utherford 
ijl99,^ for details about this process. iLeake et al.l (l2013[ l 
has shown that under chromospheric conditions this pro¬ 
cess increases the collisional coupling between ions and 
neutrals. It is therefore logical to think that this mecha¬ 
nism can be also relevant under prominence conditions. 

In the present Le tter we follow an appr oach similar to 
the one described in iTerradas et al.l (j2013[ l for a fully ion¬ 
ized plasma and is based on solving the time-dependent 
problem in a two-dimensional geometry with a magnetic 
held that incorporates dips. The main difference is in the 
injection of a plasma dominated by neutrals at the core 
of the prominence. With the two-fluid approximation 
the effects of ion-neutral collisions and charge exchange 
collisions are properly incorporated to the model. These 
mechanisms provide a redistribution of momentum of the 
species and change the force balance in the system. This 
is the crucial step to accomplish a new stationary state 
that is in dynamical equilibrium. 

2. BASIC TWO-FLUID EQUATIONS 

We use the most simplified version of the two-fluid 
equations including ion-neutral collisions and charge ex- 
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change collisions. For the ion-electron fluid we have the 
equations of continuity, momentum, pressure, and mag¬ 
netic induction. 
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We have used the simplest version of the Ohm’s law in 
which Ohmic resistivity, the Hall and battery terms are 
ignored since they are very small under typical promi¬ 
nence conditions. 

For the fluid composed of neutrals we have 


^ = -V-(p„v„), (5) 
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9pn , _ 
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In these equations all the variables have the usual mean¬ 
ing. We define the drift velocity as vd = Vi — Vn. The 
coupling between the two fluids is through the param¬ 
eters ain and A'j?- Tti® fi^st parameter represents the 
friction co efficient due to co l lisions between ions and neu - 
trals fe.g.. lBraginskiil[l965l : iChaoman fc Cowlind 1197011 . 
For hydrogen and assuming that ions and electrons have 
the same temperature the friction coefficient is given by 
the following expression 
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where rrip is the proton mass and Cin is the momen¬ 
tum transfer cross section, taken to be equ al to 10“^®m^ 
in this work (see iVranies fc Krstid 120131 1. The ther¬ 
mal speed of the different species is given by ut = 
y/2kBT/rup where T is the temperature. It is useful to 
define the collision frequency between neutrals and ions 
as z/ni = ttin/Pn- The maximum collision frequency in 
our configuration is Vni = 118 Hz. 

The second parameter, is responsible for the 

charge exchan ge collisions and is given by th e following 
expression (see iPauls et 1111119951: lMeieil[20i^ 
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For interactions between protons and neutral hydro- 
gen the cross section, CTcx, is around 10“^®m^ (see 
iMeier fc Shumlaljl2012[l . Hence, the cross sections asso¬ 
ciated to elastic collisions and charge exchange collisions 
are approximately the same. 

The terms Win, Wni, Wj^, and W“ in the pressure 
equations contain frictional heating and thermal transfer 
between the two populations associated to the two fric¬ 
tional mechanisms. The e xplicit fo r m of W n = —Wni can 
be fou nd in, for example, iDraind ()1986[ 1 or iLeake et al.l 
(j2014ll . The expressions for Wj^ = — W“, which are 
complex because the charge exchange ph ysics i s more 
com plicated, can be fo und in the works of iMeieil (1201111 
and iLeake et al.l (|2012ll . Because of the thermal transfer 
between the species included in these terms the temper¬ 
atures of ions and neutrals will tend to be equal on short 
time scales. 

Equations (|8])-(l9l) indicate that the friction coefficients 
depend essentially on the density and pressure of the 
individual species (i.e., the thermal speed). Since the 
coefficients depend on the local plasma parameters they 
also vary in space and time during the evolution, allowing 
us to model self-consistently the coupling between ions 
and neutrals. 

From Eqs. and ([n|) we have that the force between 
the fluids is proportional to the difference in velocity be¬ 
tween ions and neutrals, vd. Static stationary solutions 
are simply not possible. If neutrals are sustained against 
gravity it means that there must be a flow that provides 
the required frictional restoring force. The question is 
whether the flow produces an important draining effect 
of neutrals from the prominence body. 

3. INITIAL CONDITIONS AND METHOD 

In our model the density of ions and neutrals is con¬ 
structed using a static background plus a localized en¬ 
hancement representing the prominence. For the back¬ 
ground we assume a stratified atmosphere for ions with 
a scale height characteristic of 10® K, with gravity {g = 
274 ms“^) pointing in the negative z—direction. 

The prominence is represented by a la rge density en¬ 
hanc ement respect to the background ("see lTerradas et ahl 
IWl . imposing that its maximum is 100 times larger 
than the coronal density. This enhancement is assumed 
to be composed by 75% of neutrals and 25% of ions. 
These per centages are consistent with t he ionization 
models of iGouttebroze fc Labrossel (1200911 (see middle 
panel of their Figure 1). Regarding the spatial struc¬ 
ture of the enhancements for both ions and neutrals is 
given by a two-dimensional Gaussian with a characteris¬ 
tic width a and b in the x and z—directions, respectively. 
We use that a = 0.3 Lp and b — 0.6 Lp. 

For the magnetic configuration we use a force- 
free q uadrupolar structure described in iTerradas et al.l 
(j2013ll . Magnetic dips are present in this topology 
(see Figure [T]) and are important to achieve a sustained 
prominence. A small shear component is included in the 
model to avoid the plasma-/? to be infinite at the base 
of the corona. We use a magnetic field of 50 G at this 
level which gives 10 G at the core of the prominence (the 
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plasma-/3 is around 0.1) and the changes in the magnetic 
held during the evolution are rather small. 

The system of time-dependent nonlinear equations 
given by Eqs. 0-0 together with the initial conditions 
are solved in two-dimensions. We choose the following 
normalizing values for the length, Lq = 10^ m, and den¬ 
sity, po = 5.2 X 10“^^kgm“^. The normalizing value for 
velocity is Cso = 1-66 x 10®ms“^. The simulation domain 
extends from —6Lq to 6Lq in the a;—direction and from 0 
to 8Lq in the z—direction. We use 200 x 200 grid points. 
Line-tying conditions are imposed at z = 0 representing 
the base of the corona while extrapolated conditions are 
enforced at the lateral and upper boundaries. The inclu¬ 
sion of the frictional terms imposes a strong limitation 
on the time steps required to properly resolve the short 
diffusion time scales. 

4. RESULTS 

A representative simulation with the parameters given 
in the previous section is described here. For a better 
understanding of the results we first neglect charge ex¬ 
change collisions, which will be included later, and con¬ 
centrate on elastic scattering between ions and neutrals. 

4.1. Ion-neutral collisions 

From the analysis of the time-dependent simulations 
we find several evidences about the support of neutrals 
located at the core of the prominence. As an example, we 
have plotted in Figure [1] top panel, the two-dimensional 
distribution of the two plasma densities after only 4.8 min 
of evolution. The ion-neutral fluid and neutrals are es¬ 
sentially superimposed, and the compact spatial shape of 
the densities found in this figure is very similar to that at 
the beginning of the simulation meaning that the mass 
redistribution is rather small in the configuration. On 
the contrary, in Figured] bottom panel, the results of the 
same simulation are shown but with a reduced maximum 
value for ain (ten times smaller than in top panel). The 
behavior of the system is very different with respect to 
the previous case. Neutrals are falling quite quickly and 
diffusing across the magnetic field. A sustained situation 
for neutrals is not achieved in this case. 

We return to the case studied in Figuredjtop panel. We 
have allowed the system to evolve for more than 10 hours, 
and the geometrical changes found at the prominence 
core are very small. In fact neutrals are moving down¬ 
wards but with such small velocity that makes, to all 
practical purposes, difficult to distinguish from a static 
situation. These results clearly indicate that a sustained 
partially ionized plasma is a feasible configuration. The 
crucial point is that the friction coefficient is very large 
at the core of the prominence and the coupling between 
ions and neutrals is very strong. Essentially the multi- 
component plasma behaves as a single fluid at the promi¬ 
nence body and this means that the velocities for the 
two plasma components are very similar. Thus, the drift 
velocity is rather small, but since the frictional coeffi¬ 
cient is very large the frictional force, which is simply 
the product ainVD, is sufficient to balance gravity. This 
drift velocity, vd, is always positive and therefore the 
frictional force counterbalances the gravity force point¬ 
ing downwards and acting on neutrals. The gas pressure 
at the center of the prominence has a maximum, and 
the pressure gradient is zero. The first term on the right 
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Fig. 1. — Magnetic field (blue lines) and density of the ion- 
neutral fluid (blue-green colors) and neutrals (yellow-pink col¬ 
ors) at a given time {t = 4.8min). In the top panel the 
simulation corresponds to a maximum collision frequency of 
Ujii = 8.5 Hz, approaching to the strong coupled case. For the 
bottom panel i/ni = 0.8 Hz, the weak coupling between the 
ionized plasma and neutrals produces the diffusion of neu¬ 
trals across the magnetic field. 

hand side of Eq. is therefore negligible since p-aV^z 
quite small (we concentrate only on the z—component). 
The remaining terms provide the force balance in the 
z—direction, mean ing that k, q/v r, ^ in agr e ement 
with the results of iGilbert et al.l (j2002l l: GilbertI (l2011f l 
based on simple assumptions and the balance in the mo¬ 
mentum equation. This means that the drift velocity is 
inversely proportional to the neutral-ion frequency. So if 
the system is in a stationary state under balance of forces 
the stronger the coupling between the two-fluids (higher 
neutral-ion frequency) the smaller the vertical compo¬ 
nent of the drift velocity. In fact, the ionized plasma is in 
this regime essentially static (viz ~ 0), meaning that the 
downward velocity of neutrals is simply Vnz ~ —gl^ni- 
Now we turn out attention to the ionized plasma and 
consider Eq. ©• For the same reasons as before the gra¬ 
dient of the inertial and pressure terms is negligible but 
the magnetic term is very relevant. It provides the bal¬ 
ance against the gravity term and the frictional force that 
now is pointing downwards. To have force balance the 
magnetic force must be equal to Pig-\-ain udz, which using 
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the previously estimated value for udz, it reduces to the 
simple and illustrative expression g (pi + pn). The phys¬ 
ical interpretation is clear. The ionized fluid is always 
supported by the magnetic field but if a neutral compo¬ 
nent is included in the model the only way to have a new 
balance of forces is to increase the restoring magnetic 
force in such a way that the weight of the two plasma 
components is compensated by the deformation of the 
magnetic field. This result is only valid when the two- 
fluids are strongly coupled, and this is true at the core 
of the prominence body. 

We have explored the dependence of the results on the 
crucial magnitude of our model, namely Oin or equiva¬ 
lently Vni- We have artificially imposed the maximum 
of this parameter and have performed different numer¬ 
ical experiments. The results are shown in Figure [2l 
When i/ni is around three orders of magnitude below the 
maximum value given by Eq. m (z/ni = 118Hz), neu¬ 
trals are weakly coupled with the ion-electron fluid and 
they essentially fall down rapidly toward the base of the 
corona (see Figure [T] bottom panel). However, when this 
parameter is raised a new dynamical equilibrium is ac¬ 
complished with a reduced drift velocity. We have been 
able to check using the numerical results that the ex¬ 
pression Vnz ~ —givni IS Satisfied (compare circles with 
the dashed line), which is a further confirmation that 
the system is under a dynamical balance. In fact, due 
to numerical issues related to the size of the time-step in 
the simulations when the collision frequency is large, we 
have been forced to fix a maximum collision frequency 
below the maximum value inferred from Eq. ®. Never¬ 
theless from Eigure[2]we see that we are approaching to 
the analytic results. 

Using the maximum collision frequency (i^ni = 118 Hz) 
in our configuration, and the expression for the velocity 
drift, we obtain that the draining of neutrals takes place 
at a speed of only Vmo ~ —2.3ms“^ (see the dotted ver¬ 
tical line in EigureH]). This is a very low vel ocity and is in 
good agreement with the value inferred by [Gilbert et al.l 
()2002fl (unzo ~ —3.7ms“^). Hence, the draining of neu¬ 
trals from the prominence core proceeds on very large 
time-scales which are of the order of months using the 
simple time estimate 2To/knzo| (the time required for 
the neutrals located at a height 2 Lq to reach the base of 
the corona moving at a constant velocity |unzo|)- 

In the two-fluid approach the heating rate due to the 
conversion of kinetic energy is simply ainVo^- Using the 
expression for the drift velocity we obtain PnP^ jvni. This 
magnitude at the core of the prominence (using the max¬ 
imum frequency of 118 Hz) gives a heating rate which is 
typically of the order of 10“® J m“^ s“^. Hence, heating 
by frictional collisions is quite irrelevant. 

4.2. Charge exchange collisions 

Now the terms associated to charge exchange inter¬ 
actions are activated in the simulations. These terms 
increase the frictional coupling between the species, and 
this can be easily inferred from the comparison of Eqs. m 
and ®. Neglecting the drift velocity in front of the ther¬ 
mal velocities, and assuming that the thermal velocities 
of the ion-electron and neutral fluids are the same, it 
is straight forward to find that ~ 2.72 (since the 
two processes have the same cross section). Thus, charge 
exchange collisions produce almost three times more fric- 
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Fig. 2.— Drift velocity for different values of oin (z^ni) in¬ 
ferred from the simulations at the point x = Q, z — 2Lo. 
Circles correspond to elastic ion-neutral collisions while dia¬ 
monds to elastic plus charge exchange collisions. The dashed 
line represent to the analytic approximation vuz ~ gjv^i, 
while the dot-dashed line corresponds to the approximation 
vuz ss y/(3.72z/ni). The vertical dotted line is the maximum 
collision frequency. 

tion than that of ion-neutral collisions. This means that 
the downward velocity of neutrals is further reduced by 
charge exchange interactions and the final expression is 
Vnz - -gpn/ (ctin + /?“) ~ “ff / (3-72 . To test the as¬ 
sumptions made to derive this expression we have again 
changed the maximum collision frequency and have cal¬ 
culated from the simulations the numerical value of v^z 
when both ion-neutral and charge exchange collisions are 
present. The results are overplotted in Figure [5] with di¬ 
amonds. Once more the agreement between the simu¬ 
lations and the analytic approximation (plotted with a 
dot-dashed line) is remarkable. 

5. DISCUSSION 

In spite of the simplicity of our model it contains the 
very basic ingredients to demonstrate that the support 
of neutrals in solar prominences through frictional cou¬ 
pling is a viable mechanism. We have not tried to model 
the formation process and have concentrated mainly on 
the issue of the support of neutrals. We have found that 
as long as the friction coefficient is large, and the esti¬ 
mations of this coefficient under prominence conditions 
point in this direction, it is relatively easy to find a “dy¬ 
namical equilibrium”. We have shown for the first time 
by solving the time-dependent problem that through the 
frictional coupling, magnetic forces rearrange to balance 
the gravity force acting on the joint mass of neutrals 
and ions. Interestingly, for hydrogen the friction coeffi¬ 
cient for change exchange collisions is three times larger 
than the friction coefficient for ion-neutral collisions. The 
nonstatic equilibrium is characterized by a rather small 
downflow velocity for neutrals, around 2.3 ms“^, for ion- 
neutral collisions. But when charge exchange collisions 
are added to ion-neutral collisions this downflow velocity 
is reduced up to a factor four, since the net friction is 
increased. Therefore, the draining of neutrals from the 
prominence is not an important issue when the frictional 
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coupling between the species is strong. To all practical 
purposes, it would be difficult to distinguish a static sit¬ 
uation from the steady-state situation obtained in this 
work with such low values of the drift velocity. 

Many effects have been ignored in the present work. In 
particular, the ionization fraction does not change consis¬ 
tently because photoionization and recombination have 
been neglected. It is relatively simple to include recom¬ 
bination in the two-fluid equations but photoionization 
is much more difficult. To properly address this problem 
the full NLTE (non-local thermodynamic equilibrium) 
radiative transport equations should be coupled with the 
two-fluid equations and this is out of the scope of this 
work. Viscosity and other non-ideal effects have been 
neglected. It has also been assumed that the prominence 
is an isolated system, more representative of quiescent 
prominences, ignoring dynamical effects such as strong 
flows along the magnetic field that are often observed in 
active region prominences which could modify the drain¬ 
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ing rate of neutrals. 

Finally, instead of 2D, 3D prominences models should 
be studied, since are able to develop several kinds of in¬ 
stabilities such as the magnetic Rayleigh-Taylor insta¬ 
bility (see [HillleL^ES] [Mill ElISSl^IiD[2Qill)- The 

effect of neutrals on this instability needs to be investi¬ 
gated further since the growth-rates can be significan tly 
modified (see lDiaz et ali[2012HKhomenko et al.ll201^ . 
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